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Abstract 

A search for long-lived particles is performed using a data sample of 4.7 fb _1 from proton-proton 
collisions at a centre-of-mass energy = 7 TeV collected by the ATLAS detector at the LHC. No 
excess is observed above the estimated background and lower limits, at 95% confidence level, are 
set on the mass of the long-lived particles in different scenarios, based on their possible interactions in 
the inner detector, the calorimeters and the muon spectrometer. Long-lived staus in gauge-mediated 
SUSY-breaking models are excluded up to a mass of 300 GeV for tan/3 = 5-20. Directly produced long- 
lived sleptons are excluded up to a mass of 278 GeV. tf-hadrons, composites of gluino (stop, sbottom) 
and light quarks, are excluded up to a mass of 985 GeV (683 GeV, 612 GeV) when using a generic 
interaction model. Additionally two sets of limits on fl-hadrons are obtained that are less sensitive to 
the interaction model for 7?-hadrons. One set of limits is obtained using only the inner detector and 
calorimeter observables, and a second set of limits is obtained based on the inner detector alone. 
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1. Introduction 

Heavy long-lived particles (LLPs) are pre- 
dicted in a range of theories which ex- 
tend the Standard Model (SM). Supersym- 
metry (SUSY) [Ml models allow long-lived 
charged sleptons (£), squarks (q) and gluinos 
(g). Heavy LLPs produced at the Large 
Hadron Collider (LHC) could travel with 
speed measurably lower than the speed of 
light. These particles can be identified and 
their mass, m, determined from their speed, 
fi, and momentum, p, using the relation m = 
p/y/3, with y being the relativistic Lorentz 
factor. Four different searches are presented 
in this Letter, using time-of-flight to measure 
ft and specific ionisation energy loss, dE/dx, 
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to measure /3y. The searches are optimised 
for the different experimental signatures of 
sleptons and composite colourless states of a 
squark or gluino together with SM quarks and 
gluons, called 7?-hadrons. 

Long-lived charged sleptons would inter- 
act like muons, releasing energy by ionisa- 
tion as they pass through the ATLAS detec- 
tor. A search for long-lived sleptons iden- 
tified in both the inner detector (ID) and in 
the muon spectrometer (MS) is therefore per- 
formed ("slepton search"). The results are in- 
terpreted in the framework of gauge-mediated 
SUSY breaking (GMSB) EMI with the 
light stau (f ) as the LLP. In these models a sub- 
stantial fraction of the events would contain 
two LLP candidates, and this feature is also 
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utilised in discriminating signal from back- 
ground. Direct pair production of sleptons is 
also used to interpret the data independently 
of the mass spectrum of the other SUSY parti- 
cles. 

Coloured LLPs (q and g) would hadronise 
forming 7?-hadrons, bound states composed of 
the LLP and light SM quarks or gluons. They 
may emerge as charged or neutral states from 
the pp collision and be converted to a state 
with a different charge by interactions with 
the detector material, and thus arrives as neu- 
tral, charged or doubly charged particles in the 
muon spectrometer. 

In ATLAS LLPs can be identified via the 
timing information in the muon spectrome- 
ter or calorimeters and via the measurement 
of the energy loss in the silicon pixel detec- 
tor. All these techniques are combined in 
this analysis to achieve optimal sensitivity for 
the "full-detector 7?-hadron search". In addi- 
tion searches based on only the calorimeter 
and pixel detector information ("MS-agnostic 
7?-hadron search"), and based solely on the 
pixel detector ("ID-only 7?-hadron search") are 
performed. The latter two cases are motivated 
by the limited understanding of 7?-hadron in- 
teractions in matter, in particular the possibil- 
ity that i?-hadrons are electrically neutral in the 
MS. Furthermore, these searches are sensitive 
to scenarios in which the i?-hadrons decay be- 
fore reaching the MS. In all searches the sig- 
nal particles are assumed to be stable within 
the ATLAS detector, at least to the point it hits 
the last relevant component of the subdetector 
used for detecting it. 

Previous collider searches for LLPs have 
been performed at LEP IfTTTO. HERA (2D, 
the Tevatron IMfl. and the LHC ||29M35l 



2. Data and simulated samples 

The work presented in this Letter is based 
on 4.7 fb 1 of pp collision data collected at a 
centre-of-mass energy yfs = 7 TeV in 2011. 
The events are selected online by muon trig- 
gers for the slepton search and by missing 
transverse momentum and muon triggers for 
the 7?-hadron searches. Data and Monte Carlo 
Z — » p.fi samples are used for timing resolution 
studies. Monte Carlo signal samples are used 
to study the expected signal behaviour and to 
set limits. 

The GMSB samples are generated with 
the following model parameters: number 
of super-multiplets in the messenger sector, 
N 5 = 3, messenger mass scale, m messeng er = 
250 TeV, sign of the Higgsino mass param- 
eter, sign(/i) = 1, and C grav , the scale factor 
for the gravitino mass which determines the 
f lifetime was set to 5000 to ensure that the 
f does not decay in the detector. The ratio 
of the vacuum expectation values of the two 
Higgs doublets, tan/?, is varied between 5 and 
40 and the SUSY-breaking mass scale A is 
varied from 50 to 150 TeV, corresponding to 
light f masses varying from 122.2 to 465 GeV. 
The mass spectra of the GMSB models are ob- 
tained from the Spice program ll36ll and the 
events are generated using Herwig Il3~7ll . 

The 7?-hadron samples are generated with 
gluino (squark) masses from 300-1500 GeV 
(200-1000 GeV). The pair production of 
gluinos and squarks is simulated in Pythia 
(381, incorporating specialised hadronisation 
routines [|39l - l4T1l to produce final states con- 
taining i?-hadrons. A 10% gluino-ball frac- 
tion is assumed in the gluino sample pro- 
duction. The simulation of i?-hadron inter- 
actions with matter is handled by dedicated 
Geant4 11421 1431 routines based on a generic 
model Il4"4l . All Monte Carlo events pass the 
full ATLAS detector simulation 021 S3l and 



are reconstructed with the same programs as 
the data. All signal Monte Carlo samples are 
normalised to the integrated luminosity of the 
data. 

3. The ATLAS detector 

The ATLAS detector MB is a multipur- 
pose particle physics detector with a forward- 
backward symmetric cylindrical geometry and 
near An coverage in solid angleQ The ID 
consists of a silicon pixel detector, a silicon 
micro-strip detector, and a transition radiation 
tracker. The ID is surrounded by a thin super- 
conducting solenoid providing a 2 T magnetic 
field, and by high-granularity liquid-argon 
sampling electromagnetic calorimeters (LAr). 
An iron/scintillator-tile calorimeter provides 
coverage for hadrons in the central rapidity 
range. The end-cap and forward regions are in- 
strumented with liquid-argon calorimeters for 
both electromagnetic and hadronic measure- 
ments. The MS surrounds the calorimeters 
and consists of three large superconducting 
air-core toroids each with eight coils, a system 
of precision tracking chambers, and detectors 
for triggering. 

The ATLAS trigger system is designed to 
select the events of most interest with a data- 
taking rate of about 400 Hz from a beam bunch 
crossing rate as high as 40 MHz. The first- 
level trigger (level- 1) selection is carried out 
by custom hardware and identifies detector re- 
gions and the bunch crossing for which a trig- 
ger element is found. The high-level trigger 

'ATLAS uses a right-handed coordinate system with 
its origin at the nominal interaction point in the centre 
of the detector and the z-axis coinciding with the axis 
of the beam pipe. The x-axis points from the interaction 
point to the centre of the LHC ring, and the y-axis points 
upward. Cylindrical coordinates (r, <p) are used in the 
transverse plane, <p being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the 
polar angle 8 as 77 = - In tan(#/2). 



is performed by dedicated software, seeded by 
data acquired from the bunch crossing and re- 
gions found at level- 1. The components of 
particular importance to this analysis are de- 
scribed in more detail below. 

3.1. The pixel detector 

As the innermost detector system in AT- 
LAS, the silicon pixel detector provides at 
least three precision measurements for each 
track in the region \-q\ < 2.5 at radial distances 
from the LHC beam line r < 15 cm. The 
sensors in the pixel barrel (covering the cen- 
tral |?7|-region) are placed on three concentric 
cylinders around the beam-line, whereas sen- 
sors in the end-cap (covering the high-|?7| re- 
gion) are located on three disks perpendicular 
to the beam axis on each side of the barrel. In 
the barrel (end-cap) the intrinsic accuracy is 
10 /im in the r0-plane and 115 fim in the z(r)- 
direction. The data are only read out if the sig- 
nal is larger than a set threshold. The time for 
which the signal exceeds that threshold, ToT, 
is recorded. The larger the initial signal is the 
longer this time. 

3.1.1. Pixel detector specific ionisation 
(dE/dx) measurement 
The relation between the ToT and the charge 
deposition in each pixel is measured in dedi- 
cated calibration scans and shows a good lin- 
earity. Therefore the ToT measurement is well 
correlated with the energy loss of a charged 
particle in the pixel detector. The maximum 
ToT value corresponds to 8.5 times the av- 
erage charge released by a minimum ionis- 
ing particle (MIP) for a track perpendicular to 
the silicon detectors and leaving all its ionisa- 
tion charge on a single pixel. If this value is 
exceeded, the ToT (and therefore the charge) 
is not correctly measured. In LHC collisions 
the charge generated by one track crossing the 
pixel detector is rarely contained in just one 



pixel. Neighbouring pixels are joined together 
to form clusters and the charge of a cluster is 
calculated by summing up the charges of all 
pixels after calibration correction. The specific 
energy loss dE/dx is defined as the average of 
all individual cluster charge measurements for 
the clusters associated with the track. To re- 
duce the Landau tails, the average is evaluated 
after having removed the cluster with the high- 
est charge (the two clusters with the highest 
charge are removed for tracks having five or 
more clusters). 

3.1.2. Mass measurement with the pixel detec- 
tor 

The masses of slow charged particles can 
be measured using solely the ID information 
by fitting each dE/dx and momentum mea- 
surement to an empirical Bethe-Bloch func- 
tion and deducing their fiy value. The mea- 
surable f3y range lies between 0.2 and 1.5, the 
lower bound being defined by the overflow in 
the ToT spectrum, and the upper bound by the 
overlapping distributions in the relativistic rise 
branch of the curve. This particle identifica- 
tion method [47 1 uses a five-parameter func- 
tion to describe how the most probable value 
of the specific energy loss (Mae) depends on 

fir- 

MvLifiy) = |i ln(l + (p 2 /3yD - p, (1) 

Ax MP3 

Figure [TJ (left) shows how this function 
describes data for low momentum tracks. 
Figure [TJ (right) shows the simulated pixel 
dE/dx spectra for singly-charged hypothet- 
ical 7?-hadrons of masses 100, 300, 500 
and 700 GeV. As expected, these distribu- 
tions extend into the high pixel dE/dx re- 
gion even for high momentum tracks. The 
most probable value of dE/dx for MIPs is 
about 1.2 MeVg _1 cm 2 with a spread of about 



0.2 MeVg _1 cm 2 and a slight ?] dependence, in- 
creasing by about 10% from low-77 to high-77 
regions. 

For all tracks having a reconstructed mo- 
mentum p and a measured specific energy loss 
dE/dx above the value for MIPs, a mass es- 
timate mpy = p/fiy is obtained by inverting 
Equation ([TJ). The procedure is continuously 
monitored through precise (< 1%) measure- 
ments of the mass of known particles (kaons 
and protons). For LLPs the expected dE/dx 
values are much larger than those of SM par- 
ticles, allowing for identifying them based on 
this information. 

3.2. Calorimeters 

Liquid argon is used as the active detec- 
tor medium in the electromagnetic (EM) bar- 
rel and end-cap calorimeters, as well as in the 
hadronic end-cap (HEC) calorimeter. All are 
sampling calorimeters, using lead plates for 
the EM calorimeters and copper plates for the 
HEC calorimeter. The barrel EM calorimeter 
covers the region \tj\ < 1.475 and consists of 
three layers and a pre- sampler. The EM end- 
cap calorimeter consists of three layers in the 
region 1.375 < \rj\ < 2.5 (two for 2.5 < \r/\ < 
3.2) and a pre- sampler for 1.5 < ]?7| < 1.8. The 
four layers of the HEC calorimeter cover the 
range 1.5 < \rj\ < 3.2. 

The ATLAS tile calorimeter is a cylindrical 
hadronic sampling calorimeter. It uses steel 
as the absorber material and plastic scintilla- 
tors as the active layers. It covers radii from 
2280 to 4230 mm and the 77 coverage extends 
to I77I < 1.7. The calorimeter is subdivided into 
a central barrel covering \tj\ < 1.0 and an ex- 
tended barrel covering 0.8 < \tj\ < 1.7. Both 
barrel parts are divided into 64 modules and 
the cells in each module are divided into three 
layers. 
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Figure 1: Left: Distribution of dE/dx versus charge times momentum for minimum bias collisions in a data sample 
from 2010. With tracks reconstructed for p T > 100 MeV this sample is more suitable for the calibration of the 
Bethe-Bloch function than 2011 data in which tracks had a cutoff of pj = 400 MeV. The distribution of the most 
probable value for pions (solid), kaons (dashed) and protons (dotted) are superimposed. Right: Simulated distribution 
of specific energy loss versus momentum for singly-charged hypothetical /?-hadrons of various masses. 



3.2.1. Calorimeter fi measurement 

The ATLAS tile and LAr calorimeters have 
sufficiently good timing resolutions to distin- 
guish highly relativistic SM particles from the 
slower moving LLPs. The time resolution 
depends on the energy deposited in the cell 
and also the layer type and thickness, but typ- 
ical resolutions are 2 ns for an energy de- 
posit of 1 GeV and generally better for the tile 
calorimeter. 

To ensure the highest possible timing accu- 
racy, it is necessary to calibrate the data us- 
ing particles with known speed. This calibra- 
tion applies a common shift for each run, and 
is then performed as a function of calorime- 
ter layer and cell energy. The reliability of 
such a calibration for this analysis depends 
on the assumption that the particles used for 
calibration have similar characteristics to the 
LLPs in question when depositing energy in 
the calorimeters. The analysis uses muons for 
this purpose and it is cross-checked that jets 
give a consistent result. The effect of a possi- 



ble bias in the measured time for late-arriving 
signals was tested by applying an arrival-time- 
dependent resolution function and found neg- 
ligible. The fi measurements from all cells as- 
signed to the extrapolated track are combined 
in a cell-energy-weighted average, typically 
using timing measurements from three or four 
calorimeter cells. The resolution of the result- 
ing /? measurement can be seen in Figure [2] 
(left). For data (MC) the mean J3- value is 0.983 
(0.986) and the resolution is 0.090 (0.092). 

3.3. The muon detectors 

The MS forms the outer part of the ATLAS 
detector and detects charged particles exiting 
the calorimeters and measures their momenta 
in the pseudorapidity range \n\ < 2.7. It is also 
designed to trigger on these particles in the re- 
gion \tj\ < 2.4. In the barrel the chambers are 
arranged in three concentric cylindrical shells 
around the beam axis, while in the two end- 
caps the muon chambers are arranged in three 
wheels that are perpendicular to the beam axis. 
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The precision momentum measurement is 
performed by monitored drift tube (MDT) 
chambers, using the rj coordinate. These 
chambers consist of three to eight layers of 
drift tubes. In the forward region (2.0 < \rj\ < 
2.7), cathode strip chambers are used in the in- 
nermost tracking layer. Resistive plate cham- 
bers (RPC) in the barrel region (|//| < 1 .05) and 
thin gap chambers in the end-cap (1.05 < \rj\ < 
2.4) provide a fast level- 1 trigger and measure 
both the rj and (p coordinates of the track. 

3.3.1. fi measurement in the MS 

The default reconstruction of particles in the 
MDT chambers [48] relies on the assumption 
that they travel with the speed of light (fi = 1). 
To improve the track quality for slow LLPs, 
the individual track segments can be recon- 
structed with different values for f3. The ac- 
tual /? of the particle is estimated from the set 
of segments with the lowest^ 2 . In a succes- 
sive combined track re-fit, including ID and 
MS hits, the particle trajectory is estimated 
more accurately. The time-of-flight to each 
tube is then obtained using the difference be- 
tween the time-of-flight corresponding to the 
re-fitted track position in each tube and the 
time actually measured. By averaging the fi 
values estimated from the time-of-flight in the 
different tubes an improved MDT fi estimation 
can be achieved. 

The RPC have an intrinsic time resolution 
of ~1 ns while the digitised signal is sampled 
with a 3.12 ns granularity, allowing a measure- 
ment of the time-of-flight. In the RPC, f3 is first 
calculated separately for each hit from the in- 
dependent position and time measurements. A 
single f3 estimate is obtained by averaging the 
/3 values from all the hits. 

By definition, in a perfectly calibrated de- 
tector, any energetic muon coming from a col- 
lision at the interaction point will pass the de- 
tector at to = 0. The to distributions in the dif- 



ferent detector systems are measured and their 
means used to correct the calibration. The ob- 
served width of these distributions after cor- 
rection is used as the error on the time mea- 
surement in the f3 fit and to smear times in the 
simulated samples. The time resolution in the 
MS (about 3 ns) does not depend on the arrival 
time. 

3.4. Combining /3 measurements 

There are three possible f3 measurements 
per track, from the RPC, the MDT and the 
calorimeters. The J3 measurements from the 
different detectors are only used if f3 > 0.2 (the 
limit of the sensitivity) and if they are consis- 
tent internally, i.e. the^ 2 probability of the av- 
erage between hits is reasonable (calorimeter) 
or the RMS of the measurement is consistent 
with the expected errors (MS). Measurements 
that are accepted are combined in a weighted 
average. The weights are obtained from the 
calculated error of each measurement multi- 
plied by the pulls of the /3 distributions for 
muons from Z boson decays. 

Since J3 is estimated from the measured 
time-of-flight, for a given resolution on the 
time measurement, a slower particle has a bet- 
ter /3 resolution. Prior to the estimation, the 
timing values of the hits in the MC samples 
are smeared to reproduce the resolution mea- 
sured in the data. Figure [2] (right) shows the 
J3 distribution for selected Z — » fi/u candidate 
events in data compared to simulation after 
having smeared the hit times in simulation as 
described above. The data are well reproduced 
by the smeared MC distribution: the mean val- 
ues are 1 .000 and 0.996 and the resolution val- 
ues are 0.035 and 0.033 for data and MC, re- 
spectively. The smearing mechanism repro- 
duces the measured muon distribution. The 
same time-smearing mechanism is applied to 
the signal Monte Carlo samples. 
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Figure 2: Distribution of /? for the calorimeter (left) and combined measurements (right) obtained for selected Z — > //// 
decays in data and MC simulation. The typical resolutions are quoted in the figures. 



4. LLP candidate and event selection 

4.1. Trigger selection 

This analysis is based on events collected by 
two main trigger types: single-muon and miss- 
ing transverse momentum triggers. 

4.1.1. Single-muon trigger 

The muon trigger and its performance in 
2011 data are described in detail in Refer- 
ences [|49l[50l . This analysis uses un-prescaled 
muon triggers with a p T threshold of 18 GeV. 
Offline muons are selected with p T > 50 GeV, 
well above the trigger threshold. 

Level- 1 muon triggers are accepted and 
passed to the high-level trigger only if as- 
signed to the collision bunch crossing. Late 
triggers due to late arrival of the particles are 
thus lost. The trigger efficiency for particles 
arriving late at the muon spectrometer is dif- 
ficult to assess from data, where the major- 
ity of candidates are in-time muons. This ef- 
ficiency is obtained from simulated 7?-hadron 
and GMSB events passing the level- 1 trigger 



simulation. The muon triggers are found to be 
efficient for GMSB signatures, which contain 
two typically high-/? LLPs that reach the MS, 
and additional muons stemming from neu- 
tralino decays. The trigger efficiency increases 
roughly linearly from zero at fi = 0.62 to its 
full potential of about 90% at /3 = 0.82 for 
sleptons and 7?-hadrons that reach the MS. The 
estimated trigger efficiency for GMSB slep- 
ton events is between 70% and 85%. Muon 
triggers are less efficient for i?-hadron events, 
where one or both of the 7?-hadrons may be un- 
charged as they enter the MS and /S is typically 
low. 

4.1.2. Missing transverse momentum trigger 
Since gluinos and squarks are produced via 
the strong interaction, 7?-hadron events often 
contain high-p T jets from QCD radiation. The 
modest energy depositions of the 7?-hadrons 
themselves combined with these jets natu- 
rally give rise to missing transverse momen- 
tum (magnitude denoted as E™ 1SS ). 

The i?-hadron analyses use missing trans- 



verse momentum triggers utilising only 
calorimeter information ll50l I5T1 with thresh- 
olds between 60 and 70 GeV (a full description 
of the ATLAS trigger system is given in 061). 
For the full-detector search the single-muon 
trigger described in the previous section is 
used in addition to the E™ 1SS trigger. Unlike 
the muon trigger, there is no loss of efficiency 
for the £™ ss triggers when i?-hadrons arrive 
late. The efficiency of the £™ ss triggers de- 
creases with increasing 7?-hadron mass as the 
relative contribution of ^-initiated processes 
decreases, though it remains between 15% 
and 20% for heavy i?-hadrons. 

4.2. Offline selection 

Two different signal types are studied: slep- 
tons and i?-hadrons. Given their different ex- 
pected interactions with the ATLAS detector, 
two dedicated LLP candidate and event selec- 
tions are used as described in the following 
sections. 

4.2.1. Common event selection 

Collision events are selected by requiring 
a good primary vertex, with at least three ID 
tracks, and with requirements on the position 
of the reconstructed primary vertex. The pri- 
mary vertex is defined as the reconstructed ver- 
tex with the highest £ P 2 T of associated tracks. 
Events recorded during a time where a prob- 
lem was present in one of the relevant subde- 
tectors are rejected. Cosmic -ray background 
is rejected by removing tracks that do not pass 
close to the primary vertex in z. Candidates 
with an ID track with \z^ k - z v tx \ > 10 mm or 
\d l Q k \ > 2 mm are removed, where zj, rk is the z 
coordinate and d% k is the transverse impact pa- 
rameter at the distance of closest approach of 
the track to the primary vertex. Both require- 
ments are tightened to 1 .5 mm for the ID-only 
analysis. For the analyses involving the muon 
spectrometer, events with cosmic rays are also 



rejected by a topological cut on any two candi- 
dates with opposite rj and <f> (\q\ + 772I < 0.005 
and - fa\ -n\< 0.005). 

4.2.2. Slepton candidate selection 

The analysis searching for sleptons requires 
two muon candidates in each event, because 
two sleptons are produced, and both have a 
high probability of being observed in the MS. 
Two sets of selection criteria are applied. A 
loose selection with high efficiency is used to 
select candidates in events where there are two 
LLP candidates. In events where only one can- 
didate passes the loose selection, that candi- 
date is required to pass a tight selection. 

Candidates in the loose slepton selection are 
required to have p? > 50 GeV. The pi mea- 
surements in the ID and MS are required to be 
consistent, so that the difference between the 
ID p T and the combined fit p T does not exceed 
a half of their average. Each candidate is re- 
quired to have \rj\ < 2.5. Any two candidates 
that combine to give an invariant mass close to 
the Z boson mass (+10 GeV) are both rejected. 
Candidates are also required to have associated 
hits in at least two of the three super-layers of 
precision measurement chambers in the MS. 

The number of calorimeter cells plus MS 
hits contributing to the ft measurement must 
exceed the number of detector (sub)systems 
used by three. For signal LLPs, the /3 values 
are expected to be consistent between the in- 
dividual measurements, whereas for muons a 
low /3-value is typically due to a poor mea- 
surement in one detector component and thus 
not consistent between different components. 
Therefore the estimated fi is required to be 
consistent for measurements in the same de- 
tector system, based on the hit time resolu- 
tions, and the ft measurements from the dif- 
ferent detector systems are required to be con- 
sistent with each other. The different detec- 
tor system measurements of f3 are required to 



be pair-wise consistent at the 3<x level, and the 
combined to be consistent with the fiy esti- 
mated in the pixel detector within 3cr. Finally, 
in order to reduce the muon background, the 
combined J3 measurement is required to be less 
than 0.95. 

To pass the tight selection, a candidate is 
required in addition to have p T > 70 GeV, 
at least two separate detector systems measur- 
ing /3, the number of hits minus the number of 
detector systems participating in the measure- 
ment be at least twelve, and the consistency 
between ft estimates in different detector sys- 
tems be within 2cr. These cuts are optimised 
to give better background rejection. 

Finally, a mass cut is applied on the can- 
didate mass, nip = p/yfi, calculated from the 
candidate's measured momentum and fi. This 
cut depends on the hypothetical f mass and 
is different for different points in the GMSB 
model parameters space, determined by the ex- 
pected significance of the signal. For the two- 
candidate sample, both masses are required to 
be above the cut. 

The number of background and expected 
signal events above the mass cut in the two- 
candidate and one-candidate signal regions is 
used to search for the presence of sleptons and 
7?-hadrons. 

Typical efficiencies for signal events to pass 
all cuts including the mass cut are 20% for 
each of the one and two candidate event cat- 
egories, giving 40% efficiency in total. 

4.2.3. R-hadron candidate selection 

Since the i?-hadron contains light quarks 
and gluons in addition to the squark or gluino, 
the charge of the 7?-hadron can change fol- 
lowing nuclear interactions with the detector 
material. This possibility makes it difficult to 
rely on a single detection mechanism without 
some loss of detection efficiency, as a neutral 
state would go undetected until the next nu- 



clear interaction occurs. Some of the main 
hadronic states regarded in the models consid- 
ered are neutral, and it is thus natural to take 
an inside-out approach, starting from the ID 
track and adding discriminators from outer de- 
tector systems if a signal is seen along the ex- 
trapolated track. This is reflected in the three 
different i?-hadron searches; "full-detector", 
"MS-agnostic" (not considering the MS) and 
"ID-only" (relying solely on the ID); mak- 
ing successively fewer assumptions about the 
i?-hadron scattering model and lifetime. While 
the first two differ only in their estimate of 
J3 (the "MS-agnostic" uses only the calorime- 
ters) and the utilised triggers, the ID-only se- 
lection is generally more restrictive in order to 
reject the larger backgrounds. 

In the full-detector and MS-agnostic analy- 
ses, i?-hadron candidates are required to have 
a good quality ID track with p > 140 GeV and 
\n\ < 2.5. The track must not be within an 

rj - (p distance AR = ^(An) 2 + (A(f>) 2 = 0.3 of 
any jet with p T > 40 GeV, reconstructed using 
the anti-&, jet algorithm ll52l with distance pa- 
rameter set to 0.4. Furthermore, the track must 
not have any nearby (AR < 0.25) tracks with 
Pt > 10 GeV nor pixel hits shared with an- 
other track. Tracks with momenta greater than 
3.5 TeV are rejected as unphysical. The can- 
didate must have a good dE/dx measurement 
and a good estimate of fi. The uncertainty on 
ft is required to be less than 10% when mea- 
sured with the calorimeter only, and less than 
4% when the result of a combination. 

In the ID-only analysis, selection require- 
ments are further tightened. Vertices must 
have more than four associated tracks and 
the Ej lss trigger must be confirmed off-line 
(£™ ss > 85 GeV, including MS contributions) 
to ensure better background rejection. Candi- 
date 7?-hadron tracks must have more than two 
(six) pixel (SCT) hits, impact parameters com- 
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patible with the primary vertex, p T > 50 GeV 
and p > 100 GeV. The isolation cut is 
also more severe, rejecting events that have a 
track of p T > 1 GeV within AR = 0.25 of 
the 7?-hadron candidate. Additionally, tracks 
are discarded if their momentum resolution ex- 
ceeds 50%, or if they are identified as an elec- 
tron. 

In events with multiple 7?-hadron candi- 
dates, only one - randomly chosen - candidate 
is used. 

The final signal selection in the full-detector 
and MS-agnostic analyses, optimised for each 
mass hypothesis, is based on f3y and J3 requir- 
ing fiy < 1.5 - 2.0 and /3 < 0.8 - 0.9, with 
lower values for higher masses. A signal re- 
gion is defined in the m^-m^ plane for each 
mass point. The lower value of the mass cut is 
set such that it corresponds to being 2cr of the 
mass resolution below the nominal mass value. 
The upper value is set to 2.5 TeV. For the 
ID-only analysis the final selection requires 
the measured dE/dx value to exceed a certain 
threshold which has been set as function of 
n such that the rejection of MIPs is indepen- 
dent of 77. The selection efficiency for gluino 
7?-hadrons of 900 GeV mass is about 11% and 
7% in the full-detector and MS-agnostic analy- 
ses, respectively. It is about 6% in the ID-only 
analysis. 

An alternative 7?-hadron model, which is an 
extension of the triple-Regge model used to 
describe squark i?-hadrons [|44| has also been 
considered. In this model, the signal efficiency 
is 40% lower at 300 GeV and increasing to the 
same level at 900 GeV above which it is higher, 
compared to the model used throughout this 
analysis. 

5. Background estimation 

The background for both the slepton and 
the i?-hadron searches is mostly composed of 



high-p T muons with mis-measured /3 or large 
ionisation. The background estimation is de- 
rived from data in all cases. The background 
mass spectrum is estimated by calculating a 
mass from the p T spectrum of candidates and 
the measured J3 distribution of the background 
obtained from control samples. 

The estimation of the background mass dis- 
tributions relies on two assumptions: that 
the signal-to-background ratio before applying 
cuts on /3 is small, and that the J3 distribution 
for background candidates is due to the finite 
resolution of the measurement and is therefore 
independent of the source of the candidate and 
its momentum. Checks of the validity of these 



assumptions are discussed in Section 6.3 



The detector is divided into n regions so that 
the {3 resolution within each region is similar. 
The muon {3 probability density function (pdf) 
in each n region is the distribution of the mea- 
sured /3 of all muons in the region normalised 
to one. The sample used in producing the /3 
pdf is enlarged with respect to the main se- 
lection of the analysis by lowering the p T cut 
to 30 GeV and removing the Z veto, in order 
to increase the acceptance and reduce possible 
signal contamination. Similarly, a f3y pdf is 
constructed from tracks in a background dom- 
inated region. 

The reconstructed mass distribution of 
muons in different regions of the detector de- 
pends on both J3 and momentum distributions 
through m = p/yfi. The regions also differ 
in the muon momentum distribution; therefore 
the combination of momentum with random /3 
(see below) is done separately in each region 
and the resulting mass distributions are added 
together. 

5.1. Slepton search 

The background is determined by convolut- 
ing the /3-distribution of muons with the ex- 
pected muon momentum spectrum. The mo- 
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mentum spectrum is determined by selecting 
candidates that pass all selection requirements 



listed in Section 4.2.2 apart from the require- 
ments on /3 and m. Each candidate then gets as- 
signed a fi- value by drawing it randomly from 
the muon fi distribution. If the assigned /3- 
value is below the cut value it is used (to- 
gether with the measured momentum) to as- 
sign a mass to the candidate. The statistical un- 
certainty is reduced by using each muon can- 
didate multiple times. The distribution of mass 
values obtained this way gives the background 
estimate. 

5.2. R-hadron searches 

For the i?-hadron full-detector and MS- 
agnostic analyses, the momentum, j3y and 
jS are obtained by taking a p T - spectrum of 
charged particle tracks from the data and as- 
signing randomly f3y and J3 values to these 
tracks. The mass is then calculated from the 
measured p value and the assigned value of 
fiy or /?. The mass distributions are then nor- 
malised to data by scaling to a sideband out- 
side the signal region. 

For the ID-only analysis the choice of the 
control sample takes into account the non- 
negligible correlations between p, dE/dx and 
77. The ionisation dependence on the path 
length in the sensor is not linear ll53l . so the 
pixel dE/dx depends on 77; the ionisation also 
depends on the particle J3y via the Bethe-Bloch 
formula, and therefore on its momentum, un- 
til the Fermi plateau is reached; finally p and 
77 are not uncorrected. The distributions used 
for the random- sampling are derived from two 
categories of background events. They are ob- 
tained by applying modified selections which 
ensure that signal contamination is minimised. 

• A first sample ("low-ionisation") is used 
to generate the 77 and p distributions. 
This is selected in the same way as the 



event candidates, but without the require- 
ment on high ionisation. Instead, an up- 
per bound is placed on the dE/dx (at 
(1.8 MeVg _1 cm 2 )), ensuring orthogonal- 
ity with the signal selection. 

• A second background sample ("low- 
momentum") is used to generate the 
dE/dx templates. A background sample 
free of signal but with no upper bound 
on the dE/dx is obtained by consider- 
ing tracks that have a maximum momen- 
tum of 100 GeV. Specifically, tracks in 
the "low-momentum" background sam- 
ple satisfy all the event candidate require- 
ments except that the transverse momen- 
tum cut is looser, p T > 10 GeV, and the 
momentum p is required to be between 
40 and 100 GeV (where the Fermi plateau 
has already been reached). 

A large background sample consisting of 
two million p, 77, dE/dx triplets is randomly 
generated. The momentum is first generated 
according to a binned function based on "low- 
ionisation" events. Then the pseudorapidity is 
generated according to the 77(7?) binned func- 
tions based on "low-ionisation" events. Fi- 
nally the ionisation is generated according to 
dE/dx(n) binned functions based on "low- 
momentum" events. The normalisation of the 
generated background to the selected data is 
obtained by scaling the background to the data 
before the high dE/dx cut and in the region 
of the mass distribution where no signal is ex- 
pected (mass below 140 GeV). 

6. Systematic uncertainties 

Several possible sources of systematic ef- 
fects are studied. The resulting systematic un- 
certainties are summarised in Table [T] The er- 
rors given are those on the expected yields in 
the signal region. 
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6.1. Theoretical cross-sections 

Signal cross sections are calculated to next- 
to-leading order in the strong coupling con- 
stant, adding the resummation of soft gluon 
emission at next-to-leading-logarithmic accu- 
racy (NLO+NLL) ll54H60t The nominal cross 
section and the uncertainty are taken from 
an envelope of cross section predictions us- 
ing different PDF sets and factorisation and 
renormalisation scales, as described in Refer- 
ence llSTl . This prescription leads to a 5% rel- 
ative uncertainty on the expected signal nor- 
malisation in the slepton search, and a 15% to 
30% uncertainty for the 7?-hadron search, in- 
creasing with 7?-hadron mass. 

6.2. Expected signal 

The muon trigger efficiency is calculated us- 
ing the tag-and-probe technique on Z — » /i/i 
events as described in References 11491 [50l . 
The uncertainty on the single muon trigger ef- 
ficiency is estimated to be 1%. The reduction 
in the muon trigger efficiency due to late ar- 
rival of particles is estimated from simulation. 
The effect of the difference between data and 
MC simulation in time alignment of hits in the 
muon trigger system relative to the LHC clock 
is estimated by shifting the hit time of the high- 
est /? candidate in each simulated event by 4 ns. 
The difference in trigger efficiency when this 
change is applied is between 0.5% and 1.5% 
for the different GMSB samples, and a sys- 
tematic uncertainty of 1.5% is assigned. For 
7?-hadrons the systematic uncertainty is esti- 
mated in the same way and found to be 2%. 

The £™ ss trigger used for the i?-hadrons 
relies on the emission of jets. Therefore, 
the trigger efficiency depends on the amount 
of initial and final state radiation (ISR and 
FSR). To evaluate the associated uncertainty, 
1 TeV gluino pair-production samples are sim- 
ulated in Pythia 6.4.26 using the Perugia 201 1 
tune ||62"1| . setting the radiation level low and 



high. A simple threshold curve modelling of 
the trigger is applied to all three samples. The 
largest variation from the central sample is 
found to be 8.5%. 

The Ej lss triggers use calorimeter energy 
deposits to calculate the transverse energy, 
and are thus blind to muons, which therefore 
can be used for calibration and to study sys- 
tematic errors. To evaluate the trigger effi- 
ciency, the trigger turn-on curve is obtained 
by fitting the measured efficiency vs. E™ 1SS 
in Z —* fifi events where the Z has a high 
p T , both in data and simulation. Simulated 
events are re-weighted so the distribution of 
collisions per bunch-crossing match the run- 
ning conditions of the 2011 data. These ef- 
ficiency turn-on curves are then applied to 
the expected E™ ss spectrum from simulated 
i?-hadron events. The total uncertainty is es- 
timated from three contributions: the relative 
difference between the efficiencies obtained 
using the fitted threshold curves from Z — » 
ftfi data and simulation and the differences 
in efficiency obtained from independent ±l<x 
variations in fit parameters relative to the un- 
changed turn-on curve fit for both Z — » p/j. data 
and MC simulation. The total estimated E™ 1SS 
trigger uncertainty, including effects of a 10% 
variation of the E™ 1SS scale, is a 4.5% relative 
error on the efficiency for the signal. 

In the ID-only analysis, the effect of the of- 
fline Ej lss measurement uncertainty on the to- 
tal efficiency is evaluated by applying a scale 
factor of + 10% and a smearing of up to 20% 
to the missing transverse momentum [63J. The 
resulting uncertainty depends on the 7?-hadron 
mass and is between 7.3% (200 GeV) and 
4.5% (1500 GeV) of the central value. 

Differences in the selection efficiency 
between data and MC simulation for the 
7?-hadron full-detector and MS-agnostic 
searches (excluding the final selection on fiy, 
fi, momentum and mass selection) are evalu- 



ated using Z — > /i/j events. The overall relative 
uncertainty is found to be below 1.5%. The 
effect of the jet energy scale (JES) uncertainty 
on the requirement for signal candidates in 
the i?-hadron full-detector and MS-agnostic 
searches to be separated by AR = 0.3 from jets 
with pt > 40 GeV is found to be negligible. 

To verify that the signal efficiency is insen- 
sitive to the pile-up re-weighting, a 900 GeV 
gluino sample is divided into samples of high 
(A^ vtx > 8) and low (N vtx < 8) number of re- 
constructed primary vertices. The relative dif- 
ference in reconstruction efficiency is found to 
be negligible. 

The systematic uncertainty due to the track 
reconstruction efficiency and momentum reso- 
lution differences between data and simulation 
is estimated [64J to be 0.5% on the acceptance 
for GMSB events. For analyses relying solely 
on the ID for tracking, the resulting relative 
uncertainty on the signal efficiency is 1.3%. 

Differences in the pixel dE/dx and the sub- 
sequent measurement of By between data and 
simulation are considered by comparing the 
pixel dE/dx distribution for high-p T muons 
from Z — > pp events in data and MC simula- 
tion. This uncertainty is estimated to be 5%. 
In the ID-only analysis, this effect is evalu- 
ated by comparing the ionisation in data and 
simulation for MIPs and for particles with low 
By. Scale factors are measured in the two cases 
and found to be consistent. The uncertainty on 
the signal efficiency derived by applying these 
scale factors varies between 0.2% and 5.8% 
in the studied mass range, decreasing with the 
particle's mass. 

The signal B resolution is estimated by 
smearing the measured time of hits in the MS 
and calorimeter according to the spread ob- 
served in the time calibration. The system- 
atic uncertainty due to the smearing process 
is estimated by scaling the smearing factor up 
and down, so as to bracket the distribution ob- 



tained in data. A 4% (6%) systematic un- 
certainty is found in the one-candidate (two- 
candidate) GMSB signal region. The corre- 
sponding uncertainty for i?-hadrons is 3.6%. 
The uncertainty due to the timing calibration 
of the calorimeter hits is further tested by com- 
paring a calibration obtained from jet events to 
that obtained from muons. This results in a 
1% relative systematic uncertainty on the sig- 
nal yield. 

An uncertainty of 3.9% Il65ll66l is assigned 
to the integrated luminosity corresponding to 
this dataset. 

6.3. Background estimation 

The assumption that the B-pdf is indepen- 
dent of p is tested by constructing the /3-pdf 
separately in each ^-region using a low and a 
high-momentum sample. This results in dif- 
ferences of up to 10% on the background es- 
timate. To quantify the variability of the 6 
and momentum distributions within a region 
and its effect on the background estimation for 
sleptons, the detector is sub-divided into 25 rj 
regions instead of the eight used in the analy- 
sis and the background estimated with this di- 
vision. The resulting systematic uncertainty is 
6.5% for the two-candidate events and 3.1% 
for one-candidate events. The uncertainty on 
the background distribution due to the limited 
numbers of data events used to estimate the 
momentum distribution was calculated by di- 
viding the candidate sample randomly into two 
samples and estimating the background from 
each sub-sample separately. The resulting er- 
ror in the slepton search is negligible. A com- 
parison of the muon 8 distribution in inclusive 
muon events to that from Z — » up decays also 
exhibited negligible differences. The total un- 
certainty on the background estimate for the 
slepton search is 13% for two-candidate events 
and 1 1% for one-candidate events. 



The total uncertainty in the background es- 
timation for the full-detector i?-hadron analy- 
sis is estimated as follows. The (Sy and /3 pdfs 
are obtained by considering sidebands (lower 
momenta, for i?-hadrons 40-100 GeV). Simi- 
larly, the momentum pdf is obtained from side- 
bands in fiy and/or /?. Several variations of 
the size of the sidebands (and also inclusion of 
the signal region) result in 8-10% uncertainty 
estimates for each of the two variables. To 
test the momentum dependence of the {3 pdf, 
the number of rj bins is varied. The result- 
ing systematic uncertainty is 2%. To estimate 
the uncertainty on the background distribution 
coming from the limited statistics of the data 
samples used to estimate the momentum dis- 
tribution, the candidate sample is divided ran- 
domly into two samples and the background 
estimate is derived from each sub-sample sep- 
arately. The resulting uncertainty on the full- 
detector 7?-hadron search is 2% for low mass 
and 5% for high mass. The total uncertainty on 
the background estimate from the above tests 
is 15%. 

For the ID-only analysis, the statistical un- 
certainty in the background estimate domi- 
nates for most of the mass range, and is up to 
20%. Other effects contribute at most a few 
percent. These include variations from differ- 
ent binning choices for p, j] and dE/dx in the 
two background samples, and a contribution 
from pile-up effects. 

7. Results 

The mass distributions observed in data to- 
gether with the background estimate, its sys- 
tematic error and examples of signal are shown 
in Figure [3j for the two-candidate signal re- 
gion in the slepton search (top-left), the full- 
detector 7?-hadron search (top-right), the MS- 
agnostic i?-hadron search (bottom-left) and the 
ID-only 7?-hadron search (bottom-right). 



No indication of signal above the expected 
background is observed. Upper cross-section 
limits on new particles are set by counting the 
number of events passing a set of mass cuts 
optimised for a given mass point and model. 
For the ID-only analysis the full mass spec- 
trum of the background and the hypothetical 
signal is considered. Cross-section limits are 
obtained using the CL S prescription [|67l . Mass 
limits are derived by comparing the obtained 
cross-section limits to the lower edge of the lcr 
band around the theoretically predicted cross- 
section for each process. A Poisson likeli- 
hood function is used for the 7?-hadron limit- 
setting while a likelihood combining two Pois- 
son functions is used for the slepton search to 
benefit from the separation into one-candidate 
and two-candidate signal regions. 

The resulting production cross-section lim- 
its in the GMSB scenario as a function of 
the stau mass are presented in Figure [4] and 
compared to theoretical predictions. Long- 
lived staus in GMSB models with N 5 = 3, 
OT mesS en g er = 250 TeV and sign(/i) = 1 are ex- 
cluded at 95% confidence level (CL) up to a 
stau mass of 300 GeV for tan/? = 5 - 20, and 
up to stau masses of 285 GeV and 268 GeV 
for tan/3 = 30 and tan/? = 40, with expected 
limits of 305 GeV, 287 GeV and 270 GeV 
respectively. The lower limit on A is 99- 
110 TeV for values of tan/3 from 5 to 40. For 
A values around this limit, strong production 
of squarks and gluinos is suppressed due to 
their large masses. Directly produced sleptons 
comprise 38-64% of the GMSB cross-section, 
and depend only on the slepton mass. Us- 
ing the same analysis, directly produced long- 
lived sleptons, or sleptons decaying to long- 
lived ones, are excluded at 95% CL up to a f 
mass of 278 GeV (273 GeV expected) for mod- 
els with slepton mass splittings smaller than 
50 GeV. 

The resulting limits on i?-hadrons are shown 
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Figure 3: Observed data, background estimate and expected signal in the two-candidate signal region in the slepton 
search (top-left), full-detector /?-hadron search (top-right), MS-agnostic /?-hadron search (bottom-left) and in the ID- 
only /?-hadron search (bottom-right). For the slepton search only the lower of the two masses is plotted and for the 
full-detector /?-hadron search there is one candidate in the mp y overflow bin not shown in the figure. The dashed lines 
in the upper-right plot mark the lower edges of the signal region chosen for this /?-hadron mass point. 
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Table 1: Summary of systematic uncertainties (given in percent). Ranges indicate a mass dependence for the given 
uncertainty (low mass - high mass). 



in Figures [5] and [6] compared to the theoret- 
ical cross-sections. Gluino 7?-hadrons in a 
generic interaction model are excluded up to 
a mass of 985 GeV. Stop i?-hadrons are ex- 
cluded up to a mass of 683 GeV, and sbot- 
tom i?-hadrons are excluded up to a mass of 
612 GeV. The according expected limits are 
1015 GeV, 707 GeV and 645 GeV respectively. 
The MS-agnostic search yields a lower mass 
limit of 989 GeV for gluinos, 657 GeV for 
stops and 618 GeV for sbottoms, with expected 
limits of 988 GeV, 647 GeV and 615 GeV re- 
spectively. The MS-agnostic limits are inde- 
pendent of the fraction of 7?-hadrons that arrive 
charged in the MS. 

The limits from the ID-only search are 
shown in Figure [7] compared to the theoreti- 
cal cross-sections. Gluino 7?-hadrons withm < 
940 GeV, stop 7?-hadrons with m < 604 GeV 
and sbottom 7?-hadrons with m < 576 GeV are 
excluded at 95% CL, with expected limits of 



952 GeV, 620 GeV and 571 GeV respectively. 
The ID-only search has similar sensitivity also 
to i?-hadrons with lifetimes as short as a few 
ns. 

8. Conclusion 

Searches for long-lived massive particles 
were performed using data from pp collisions 
at = 7 TeV collected by the ATLAS de- 
tector in 2011, corresponding to an integrated 
luminosity of 4.7 fbr 1 . Several different sig- 
natures are considered. The data are found 
to match the Standard Model background ex- 
pectation in all signal regions. The exclusion 
limits placed for various models impose new 
constraints on non-SM cross- sections. Long- 
lived f's in the GMSB model considered, for 
tan/3 = 5 - 20, are excluded at 95% CL for 
masses up to 300 GeV, while directly produced 
long-lived sleptons, or sleptons decaying to 
long-lived ones, are excluded at 95% CL up to 
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Figure 4: Cross-section limits as a function of the f mass in GMSB models (top-left), limits on A and tan/3 in GMSB 
models (top-right) and the cross-section limits as a function of mass on directly produced sleptons (bottom). The 
top-left and bottom plots show observed and expected limits with ±l<x and ±2<r uncertainty bands, as well as the 
theoretical prediction for the cross-section is shown with a lcr band. In the top-right plot, the observed and expected 
limit are shown, with the lcr theoretical uncertainty and the lcr band, respectively. 
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Figure 5: Cross-section upper limits at 95% CL for gluino (left) and squark (right) /?-hadrons in the full-detector search 
are shown. Both the expected limit with ±lcr and ±2<x uncertainty bands and the observed upper limit are given. In 
addition the theoretical prediction for the production cross-section calculated at NLO+NLL and its uncertainty are 
drawn. 




Figure 6: Cross-section upper limits at 95% CL for gluino (left) and squark (right) 7?-hadrons in the MS-agnostic 
search, which is independent of the fraction of /?-hadrons that are charged in the MS. Both the expected limit with 
±l<r and ±2<x uncertainty bands and the observed upper limit are given. In addition the theoretical prediction for the 
production cross-section calculated at NLO+NLL and its uncertainty are drawn. 
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Figure 7: Cross-section upper limits at 95% CL for gluino (left) and squark (right) 7?-hadrons in the ID-only search. 
Both the expected limit with ±lcr and ±2cr uncertainty bands and the observed upper limit are given. In addition the 
theoretical prediction for the production cross-section calculated at NLO+NLL and its uncertainty are drawn. 



a f mass of 278 GeV for models with slepton 
mass splittings smaller than 50 GeV. 

Long-lived i?-hadrons containing a gluino 
(stop, sbottom) are excluded for masses up 
to 985 GeV (683 GeV, 612 GeV) at 95% CL, 
for a generic interaction model. Limits 
obtained independently of the fraction of 
.K-hadrons that arrive charged in the MS ex- 
clude masses up to 989 GeV (657 GeV, 
618 GeV). Furthermore, using only the in- 
ner detector, 7?-hadrons are excluded up to 
940 GeV (604 GeV, 576 GeV). 

These results substantially extend previous 
ATLAS limits Il3"2l [331, and are complemen- 
tary to the searches for SUSY particles which 
decay promptly. 
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Figure A. 8: Coloured version of Figure [T] (left): Distribution of dE/Ax versus charge times momentum for minimum 
bias collisions in a data sample from 2010. With tracks reconstructed for pi > 100 MeV this sample is more suitable 
for the calibration of the Bethe-Bloch function than 201 1 data in which tracks had a cutoff of pj = 400 MeV. The 
distribution of the most probable value for pions (solid), kaons (dashed) and protons (dotted) are superimposed. 
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Figure A.9: Left: Muon trigger efficiency in simulated GMSB events as a function of the highest (i candidate in the 
event. Right: Distribution of /3 for candidates with the highest /3 (black) and for all candidates in the events (red). 
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Figure A. 10: Examples of the may vs. nip distributions for gluino /?-hadrons with masses of 300, 700, 900 and 
1000 GeV, with the full-detector search. The figures show data, background and signal after /3, fiy and p cuts. The 
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Figure A. 11: Examples of the mp y vs. mp distributions for gluino /?-hadrons with masses of 300, 700, 900 and 
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Figure A. 12: Observed data, background estimate and expected signal in the full-detector /?-hadron search, for various 
mass hypotheses (300, 700, 900 and 1000 GeV). The ID histograms are projections of the 2D signal range with 
unconstrained mp y . 
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Figure A. 13: Observed data, background estimate and expected signal in the MS-agnostic /?-hadron search, for various 
mass hypotheses (300, 700, 900 and 1000 GeV). The ID histograms are projections of the 2D signal range with 
unconstrained mp y . 
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Figure A. 14: Efficiency times acceptance for directly produced slepton events in the one-candidate and two-candidate 
signal regions as a function of the stau mass. 
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Figure A. 15: Efficiency times acceptance for the tested /?-hadron mass hypotheses. Each point on the x-axis represent 
a hypothesis mass, the efficiency for a given hypothesis is based on the number of signal events in the counting region 
after cuts, divided by the number of generated events. The band represents the 11.4% systematic uncertainty on the 
efficiency. 
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Figure A. 16: Efficiency times acceptance for the tested 7?-hadron mass hypotheses. Each point on the x-axis represent 
a hypothesis mass, the efficiency for a given hypothesis is based on the number of signal events in the counting region 
after cuts, divided by the number of generated events. The band represents the systematic uncertainty on the efficiency, 
varying with mass as specified in Table [T] 
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Figure A. 17: Proton mass calculated from the momentum measured in the inner detector and the specific energy loss 
measured in the pixel detector. The data, covering the 2011 statistics considered for this analysis, are subdivided in 
periods of similar data taking conditions to illustrate the stability of the method. The red horizontal line represents the 
nominal proton mass value. 
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Figure A. 18: Plot of the dE/dx most probable values (MPVs) obtained from data, selected as in the "low-ionisation'' 



sample (see Section 5.2 1, as a function of rj overlapped with a fit obtained using a third order polynomial. 
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Table A. 2: Signal efficiency, number of expected background events, events observed in data, expected and observed 
limit for the one candidate (SR\) and two candidates (SR2) signal regions in various GMSB scenarios. 
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Table A. 3: Signal efficiency, number of expected background events, events observed in data, expected and observed 
limit for the one candidate (SRi) and two candidates (SR2) signal regions for direct slepton production. 
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Table A.4: Signal efficiency, number of expected background events, events observed in data, expected (with lcr 
uncertainties) and observed limit for various gluino /?-hadron mass hypotheses in the full-detector search. 
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Table A. 5: Signal efficiency, number of expected background events, events observed in data, expected (with lcr 
uncertainties) and observed limit for various sbottom /?-hadron mass hypotheses in the full-detector search. 
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Table A.6: Signal efficiency, number of expected background events, events observed in data, expected (with lcr 
uncertainties) and observed limit for various stop /?-hadron mass hypotheses in the full-detector search. 
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Table A. 7: Signal efficiency, number of expected background events, events observed in data, expected (with lcr 
uncertainties) and observed limit for various gluino /?-hadron mass hypotheses in the MS-agnostic search. 
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Table A. 8: Signal efficiency, number of expected background events, events observed in data, expected (with lcr 
uncertainties) and observed limit for various sbottom /?-hadron mass hypotheses in the MS-agnostic search. 
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Table A.9: Signal efficiency, number of expected background events, events observed in data, expected (with Icr 
uncertainties) and observed limit for various stop /?-hadron mass hypotheses in the MS-agnostic search. 
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